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OBJECTIVES 

To  carry  out,  fundamental  and  wide  ranging  research  investigations  involving  the  nonlinear 
wave  propagation  which  arise  in  physically  significant  systems  with  emphasis  on  nonlinear 
optics.  The  modeling  and  computational  studies  of  wave  phenomena  in  nonlinear  optics 
include  ultrashort  pulse  dynamics  in  mode-locked  lasers,  localized  modes  in  waveguide  arrays 
and  photonic  lattices,  novel  phenomena  in  quadratic  nonlinear  media  and  the  dynamics  of 
dispersive  shock  waves. 

STATUS  OF  EFFORT 

The  PPs  research  program  in  nonlinear  wave  propagation  is  broad  based  and  very  at  f  ive. 
There  have  been  a  number  of  important  research  contributions  carried  out  as  part  of  the 
effort  funded  by  the  Air  Force.  During  the  period  15  March  2006  30  November  2008.  eight 

papers  were  published  in  refereed  journals,  two  book  chapters  were  published,  two  refereed 
conference  proceedings  were  published  and  nineteen  invited  lectures  were  given.  The  key 
results  and  research  directions  are  described  below  in  the  section  on  accomplishments/ new 
findings.  Full  details  can  be  found  in  our  research  papers  which  are  also  listed  at  the  end  of 
this  report. 
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Research  investigations  carried  out  by  the  PI  and  colleagues  included  the  following.  The 
modes,  dynamics  and  properties  of  mode-locked  lasers  which  are  used  to  create  ultra-short 
pulses  were  analyzed.  Titaniumisapphire  (Tiisapphire  or  Ti:s)  lasers  are  often  used  to  pro¬ 
duce  ultrashort  pulses  on  the  order  of  a  few  femtoseconds.  There  are  other  mode-locked 
lasers  which  produce  ultrashort  pulses,  such  as  Sr-Forsterite,  fiber  lasers,  and  Chromium- 
doped  lasers.  Ti:s  lasers  are  known  to  have  outstanding  characteristics.  These  inode-locked 
lasers  can  be  used  to  generate  a  regularly  spaced  train  of  ultrashort  pulses  separated  by  one 
cavity  round-trip  time.  A  typical  mode-locked  laser  system  consists  of  a  Tiisapphire  crystal 
which  exhibits  a  nonlinear  Kerr  response  and  has  a  large  normal  group- velocity  dispersion 
(GVD).  This  requires  a  set  of  prisms  and/or  mirrors  specially  designed  to  have  large  anoma¬ 
lous  GVL)  in  order  to  compensate  for  the  normal  GVD  of  the  crystal.  Recent  experiments 
conducted  at  the  University  of  Colorado,  in  collaboration  with  our  group,  demonstrated 
that  such  lasers  can  be  approximated  by  dispersion-managed  systems  and  the  intra-cavity 
pulse  was  found  to  be  described  by  a  dispersion-managed  soliton.  Improved  mathematical 
models  reflecting  gain  and  loss  mechanisms  are  being  developed;  these  models  contain  gain 
and  filtering  terms  saturated  by  energy  and  a  loss  term  saturated  by  power.  The  new  models 
describes  t  he  mode-locking  and  dynamics  of  solitons. 

A  characteristic  of  short  pulse  lasers  is  the  carrier-envelope  phase  (CEP)  slip  which  is  the 
change  in  phase  between  carrier  and  envelope  from  pulse  to  pulse  in  the  pulse  train.  This 
is  the  phase  slip  that  the  intra-cavity  pulse  accumulates  over  one  cavity  round-trip  before 
being  emitted  from  the  output  coupler.  The  intra-cavity  slip  is  affected  by  the  nonlinearity 
and  dispersion  of  the  cavity.  Control  of  the  phase  slip  allows  researchers  to  stabilize  trains  of 
ultrashort  pulses  which  are  useful  for  applications.  Improved  understanding  of  the  phase  slip 
and  noise  characteristics  will  help  experimentalists  improve  the  characterisitcs/stabihzation 
of  the  pulse  train. 

Investigations  of  pulse  propagation  in  photonic  lattices  were  carried  out.  There  has  also 
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been  important  experimental  research  on  discrete  optical  wave-guides  and  the  propagation  of 
their  nonlinear  modes.  Experimentalists  have  been  able  to  construct  one  and  two  dimensional 
lattices  by  interfering  laser  beams.  This  all-optical  technique  is  a  significant  advantage 
over  prior  methods  in  which  the  background  lattice  structure  was  developed  by  mechanical 
fabrication. 

Early  experimental  observation  of  one-dimensional  nonlinear  lattice  modes  in  corrugated 
optical  waveguide  arrays  demonstrated  that  at  sufficiently  high  power,  a  laser  beam  could  be 
self-trapped  inside  the  waveguide.  This  demonstrated  the  formation  of  a  lattice  or  discrete 
soliton.  Importantly,  such  waveguides  can  be  constructed  on  extremely  small  scales  and,  as 
mentioned  above,  recently  have  been  constructed  by  all-opitcal  means.  In  turn,  such  non¬ 
linear  waves  in  waveguide  arrays  have  attracted  special  attention  due  to  their  realizability. 
Lattice  nonlinear  Schrodinger  equations  provide  excellent  models.  We  developed  asymp¬ 
totic  and  computational  met  hods  which  describe  observed  phenomena  and  found  localized 
pulse  solutions  to  two-dimensional  optical  lattices  with  both  regular  and  irregular  lattice 
backgrounds.  Irregular  lattice  backgrounds  include  vacancy  defects,  edge  dislocations  and 
quasi-crystal  structures. 

Research  involving  quadratic,  or  so-called  \(2L  nonlinear  optical  materials  has  led  to 
a  novel  asymptotic  systems  of  equations.  Detailed  calculations  indicate  that,  in  certain 
parameter  regimes,  there  are  stable  localized  pulse  solutions.  In  other  cases,  the  equations 
have  unstable  and  singular  solutions.  The  possibility  of  such  singular  solutions  indicates 
situations  when  extreme  damage  to  the  underlying  optical  crystal  is  possible. 

Dispersive  shock  waves  (DSW’s)  were  investigated.  With  experimentalists  at  the  Univer¬ 
sity  of  Colorado,  dispersive  shock  “blast”  waves  and  interactions  were  studied.  The  relevant 
analytical  approximations  and  theory  for  DSW’s  were  developed  and  results  were  compared 
with  experiment  and  computation;  key  differences  between  DSW  and  viscous  shock  waves 
were  also  described.  DSW’s  arise  in  nonlinear  optics  and  many  other  areas  of  physics. 


3 


ACCOMPLISHMENTS/NEW  FINDINGS 


Dynamics  of  ultra-short  laser  pulses  and  frequency  combs 

Research  breakthroughs  over  the  past  few  years  with  mode-locked  lasers,  such  as  Tiisapphire 
lasers,  have  enabled  scientists  to  generate  regularly  spaced  trains  of  ultrashort  pulses,  which 
are  separated  bv  one  cavity  round-trip  time.  Fig.  ??  below  shows  a  schematic  of  a  mode- 
locked  Ti:sapplnre  laser  and  the  emitted  pulse  train.  Typical  values  for  a  Tiisapphire  mode- 
locked  laser  are  pulse  width:  r  =  10  fs  =  10~14  sec  and  repetition  time:  Trep  =  10  11s 
=  10“8  sec. 


Figure  1:  Tiisapphire  laser  (left)  and  the  emitted  pulse  train  (right). 


Figure  2:  Schematic  of  a  pulse  train  (left)  and  its  spectrum  frequency  comb  (right ). 

Associated  with  the  spectrum  of  the  pulse  train  is  a  frequency  comb,  whose  frequencies 
are  separated  by  the  laser’s  repetition  frequency  /rep  =  =  100  MHz:  u;rep  —  2~fnv:  see 
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Fig.  ??  above.  In  this  figure  Tn , 07l  are  the  center  time  and  phase  of  the  nth  pulse.  In  the 
absence  of  noise,  the  pulse’s  spectrum  determines  the  bandwidth,  while  the  repetition-l  ime, 
rr,p  =  Tn+\  —  Tn,  and  overall  phase  slip,  A <p  =  <t>n+\  —  r/>7l,  determine  the  comb  function 
in  frequency  space.  The  frequency  of  the  fc’ th  comb  line  (enumerated  around  the  center 
frequency)  is  04  =  kwrep  +  u)Q,  where  urep  =  2n/Trep  and  uj0  =  A0o;rep/27r  are  the  repetition 
and  offset  frequencies  which  arc  depicted  in  Fig.  ??.  The  linewidth  (uJ\/2  in  the  inset)  is  the 
FWHM  (full  width  half  maximum)  of  the  comb  function  around  each  comb  frequency.  In 
the  deterministic  case,  for  a  large  number  of  pulses  N  >>  1  the  linewidth  can  be  estimated 
to  be  uj\ j2  =  0{jff — )•  Additional  noise  leads  to  jitter  in  the  center  time,  Tn?  and  phase,  on, 
which  in  turn  broadens  the  FWHM  of  the  comb  lines. 

As  discussed  above,  mode- locked  lasers  such  as  Ti:sapphire  (Ti:s)  laser  systems  (see 
Fig.  ??)  can  generate  trains  of  optical  pulses,  whose  spectrum  consists  of  frequency  comb 
lines  (see  Fig.  ??).  Important  progress  in  the  development  of  extremely  stable  optical  os¬ 
cillators  has  been  made  possible  by  the  use  of  controlled  femtosecond  frequency  combs. 
Extremely  stable  frequency  combs  have  been  generated  by  Tins  laser  systems,  but  others 
such  as  Sr:Forsterite  lasers  are  also  being  studied  intensively.  We  have  been  working  with 
faculty  in  the  Department  of  Physics  at  the  University  of  Colorado  on  this  research. 

In  our  research  we  have  been  studying  a  distributed  dispersion-rnanaged,  power  energy 
saturation  (PES)  model  system.  For  a  pulse  with  amplitude  u(z,  i),  power  P(z,l)  —  |u  2, 
and  energy  E{z)  =  J  ^  |u|2  d£,  which  is  propagating  in  the  z  direction,  our  normalized 
equation  takes  the  form 

il 
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where  the  constant  parameters  ry,  r,  /,  Es at,  Psat  are  positive.  The  first  term  011  the  right 
hand  side  represents  saturable  gain,  the  second  is  nonlinear  filtering  (r  ^  0)  and  t  he  third  is 
saturable  loss.  This  model  generalizes  the  well-known  master  laser  equation.  A11  important 


observation  is  that  when  the  loss  term  is  approximated  in  the  weakly-nonlinear  regime  by  a 
first  order  Taylor  polynomial  we  obtain  the  master  laser  equation.  Hence,  the  master  laser 
equation  is  included  in  the  power  saturated  model  as  a  first  order  approximation. 

The  dimensionless  governing  equation  is  obtained  by  transforming  to  the  following  vari¬ 
ables: 

z'  =  z/z,,t'  =  t/t,,u  =  E/ \J~P, 

where  E  is  the  envelope  of  the  electromagnetic  field  t+,I\  are  the  characteristic  time  (pro¬ 
portional  to  pulse  width)  and  power  respectively.  We  take  z*  —  l/(P*7o)f  where  70  (in  units 
1/MW-mrn)  is  the  nonlinear  coefficient  in  the  laser  crystal,  g ,r  and  /  have  all  been  scaled 
by  z*  and  the  normalized  dispersion  is  given  by 

d{z)  =  -k"{z)/kl  k"  =  /*/*.. 

where  kn  is  the  GVD  (in  units  fs2/imn).  In  normalized  units  we  find 

d{z)  =  (d)  +  A(^, 

with  (cl)  being  the  average  dispersion  (net  GVD),  A  the  deviation  from  the  average  GVD, 
and  lc  the  normalized  laser-design  map  length.  Usually  one  considers  a  two  step  dispersion 
map  where  A  j,  j  =  1,2  are  taken  to  be  constant  in  the  mirror-f  prism  components  ( j  = 
1)  and  crystal  (j  —  2).  The  map  length  over  which  the  anomalous  dispersion  occurs  is 
0lr  where  0  <  0  <  1;  typically  we  take  6  =  1/2  in  Ti:sapphire  laser  applications.  We 
also  introduce  the  map  strength  parameter  s,  which  is  proportional  to  the  area  under  the 
dispersion  map,  s  =  ^[0  A  \  —  (1  —  9)  A2].  In  addition  to  dispersion-management,  we  also  have 
nonlinear- management  in  this  laser  model  .  Here  nonlinear-management  means  that  n  —  1 
(transformed  from  70  in  dimensional  variables)  inside  the  Ti:sapphire  crystal  and  n  —  0 
outside  the  crystal;  that  is,  we  assume  linear  propagation  inside  the  prisms  and  mirrors. 
In  Fig.  ??  below  a  schematic  diagram  of  the  normalized  dispersion-  managed  configuration 
associated  with  a  mode-locked  Ti:sapphire  laser  system  is  shown. 
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Figure  3:  Figure  shows  a  typical  dispersion-managed  configuration. 


In  our  earlier  research  investigations  in  liber  optics  we  derived,  based  on  the  asymptotic 
procedure  of  multiple  scales,  a  nonlinear  integro-differential  equation  (not  given  here  due 
to  space  considerations)  which  governs  the  dynamics  of  dispersion-managed  pulse  propaga¬ 
tion.  This  governing  equation  is  referred  to  as  the  dispersion-managed  nonlinear  Schrodinger 
(DMNLS)  equation.  When  there  is  no  gain  or  loss  in  the  system,  for  strongly  dispersiou- 
managed  systems,  the  DMNLS  equation  plays  the  role  of  the  “pure"  NLS  equation  which 
is  the  relevant  averaged  equation  when  there  is  either  small  or  no  dispersion-management. 
When  gam  and  loss  are  included  as  in  the  PES  equation  we  have  a  modification  of  the  “pure" 
DMNLS  equation. 

With  or  without  dispersion-management  the  PES  equation  naturally  describes  the  locking 
and  evolution  of  pulses  in  mode-locked  lasers  that  are  operating  in  the  soliton  regime.  To 
describe  our  research  in  more  detail,  we  fix  typical  values  of  the  parameters  vary  the  gain 
parameter  g  and  the  map  strength  s.  When  g  <  g*(s)  no  localized  solution  is  obtained; 
i.c.  in  this  case  the  effect  of  loss  is  stronger  than  a  critical  gain  value  and  the  evolution  of 
a  Gaussian  profile  decays  to  the  trivial  solution.  Conversely,  when  g  >  g*(s),  there  exists 
a  single  localized  solution,  u  =  Uo(t)exp(ifiz)  where  p,  called  the  propagation  constant. 
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is  uniquely  determined  given  the  specific  values  of  the  other  parameters.  The  localized 
solutions  of  the  modified  DMNLS  equation  with  gain-loss  in  the  different  regimes  are  shown 
in  Fig.  ??;  in  this  figure  we  see  that  for  given  map  strengths  the  amplitude  of  the  pulses 
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Figure  4:  Solitons  of  the  DMNLS  equation  with  gain-loss  corresponding  to  different  map 
strengths  and  gain  parameters.  Notice  that  these  solitons  occur  only  for  a  specific  value  of 
propagation  constant,  //,  in  contrast  to  the  unperturbed  case  where  solutions  exist  for  all 
//  >  0,  thus  providing  the  mode-locking  mechanism. 

increases  with  g.  In  addition,  we  see  that  the  pulses  become  broader  as  the  map  strength, 
increases.  Blow  up  is  not  obtained,  even  under  extreme  gain.  As  mentioned  earlier,  the 
influence  of  gain-loss  on  pulses  in  this  model  is  to  create  the  mode-locking  mechanism.  The 
resulting  modes  also  correspond  closely  to  the  modes  of  the  unperturbed  equation  with  the 
same  propagation  constant  //,. 

Interestingly  it  is  known  that  only  for  a  narrow  range  of  parameters  the  master  laser 
equation  (when  the  loss  term  is  taken  to  be  the  first  two  terms  of  the  Taylor  expansion  of 
the  last,  power  saturated  term,  in  the  PES  equation)  exhibits  stable  soliton  solutions  with 
mode-locking  evolution.  Otherwise  the  solitons  are  found  to  be  unstable;  either  dispersing  to 
radiation  or  evolving  into  nonlocalized  quasi-periodic  states.  In  addition,  for  diffident  pa¬ 
rameters,  the  amplitude  can  also  grow  rapidly  under  evolution.  Thus,  the  basic  master  laser 
equation  captures  some  qualitative  aspects  of  pulse  propagation  in  a  laser  cavity;  however, 
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since  there  is  only  a  small  range  of  the  parameter  space  for  which  stable  inode- locked  soliton 
pulses  exist,  it  does  not  reflect  the  wide  ranges  of  operating  conditions  where  mode-locking 
occurs. 

The  PES  and  the  master  laser  equation  have  parameters  which  are  obtained  from  experi¬ 
ments.  When  there  is  insufficient  gam  in  Eq.  (??),  both  models  show  that  pulses  dissipate  to 
zero.  Remarkably,  a  distinguishing  feature  of  the  PES  model  is  that  even  under  large  gain, 
pulses  do  not  blow-up  nor  do  they  exhibit  instabilities.  Indeed,  when  the  gain  is  greater  than 
some  threshold  value  g  =  g*,  during  the  evolution  the  pulse  readjusts  itself  as  it  mode  locks 
into  a  stable  localized  mode,  or  soliton  solution,  which  we  refer  to  as  a  soliton  wave  attractor 
(SWA).  Furthermore,  all  pulses  that  evolve  into  SWAs  can  be  obtained  independently  using 
a  mode  finding  algorithm  which  we  developed  in  our  earlier  research  supported  by  AFOSR. 

Power  saturation  models  also  arise  in  other  problems  in  nonlinear  optics  and  are  im¬ 
portant  in  the  underlying  theory.  For  example,  power  saturation  models  are  important  in 
the  study  of  the  dynamics  of  localized  lattice  modes  (solitons,  vortices,  etc)  propagating  in 
photorefractive  nonlinear  crystals.  If  the  nonlinear  term  in  these  equations  was  simply  a 
cubic  nonlinearity,  without  saturation,  two  dimensional  fundamental  lattice  solitons  would 
be  vulnerable  to  blow  up  singularity  formation,  which  is  not  observed.  Thus  saturable  terms 
are  crucial  in  these  problems. 

Carrier-envelope  phase  slip 

As  mentioned  above,  mode-locked  lasers  can  generate  a  regularly  spaced  train  of  ultra- 
short  pulses  separated  by  one  cavity  round-trip  time.  The  phase  slip  is  the  change  of  the 
phase  offset  between  carrier  and  envelope  from  pulse  to  pulse  in  the  pulse  train  which  accu¬ 
mulates  over  one  cavity  round-trip,  before  being  emitted  from  the  output  coupler.  Fig.  ?? 
depicts  the  physical  origin  of  the  carrier  envelope  phase  (CEP)  shift.  The  intra-cavitv  slip 
is  induced  by  the  nonlinearity  and  dispersion  of  the  cavity.  The  phase  offset  <jj0  in  the 
frequency  comb  is  proportional  to  the  carrier  envelop  phase  ujq  =  A<£u;rep/27r  where  the 
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repetition  frequency  is  u;rep  =  2rr /Trep  and  the  carrier  envelop  phase  is  given  by  ^Qct:  =  A0. 


Figure  5:  The  carrier-envelope  phase,  A 4>ce  —  Ay>,  changes  during  propagation,  because 
the  envelope  propagates  at  the  group  velocity  while  the  carrier  wave  propagates  at  t  he  phase 
velocity. 

A  typical  Ti:s  laser,  such  as  the  one  depicted  in  Fig.  ??,  consists  of  a  Ti:sapphire  crystal 
that  has  a  nonlinear  Kerr  response  as  well  as  large  normal  group- velocity  dispersion  (GVD), 
and  a  set  of  prisms  and  mirrors  specially  designed  to  have  large  anomalous  GYD.  The  pump 
laser  excites  the  Ti:Sapphire  crystal,  causing  it  to  lase  and  the  pulse  undergoes  large  changes 
inside  the  Thsapphire  cavity  .  The  combined  contributions  to  the  phase  slip  depends  on  the 
nonlinear  phase  and  nonlinear  dispersion  in  the  cavity.  The  crystal  induces  a  nonlinear  effect 
as  well  as  large  normal  GVD;  the  mirror  and  prisms  induce  large  anomalous  GVD  that  nearly 
balances  the  normal  GVD  of  crystal  with  a  small  net-positive  GVD  (average  dispersion)  over 
one  round-trip.  The  pulse  bounces  between  the  mirrors  and  output  coupler  and  is  “sampled" 
every  round-trip  at  the  output  coupler  (which  transmits  only  G%  of  the  energy).  When  the 
laser  is  mode-locked  a  regularly-spaced  ultrashort  pulse  train  is  emitted  from  the  cavity.  In 
our  earlier  work  we  employed  a  perturbed  NLS  equation  without  gain/loss;  i.e.  we  took 
g  =  r  =  Z  =  0  in  Eq.  (??),  but  we  added  a  small  nonlinear  term  on  the  RHS  of  proportional 
to  n(^)(|a|2zi)f  and  a  higher  order  (third  order)  linear  dispersive  term  (see  below).  The 
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nonlinear  term  on  the  right-hand  side,  often  called  the  “shock"  term,  corresponds  to  small 
nonlinear  dispersion  arising  from  the  Kerr  effect.  The  shock  term  is  particularly  important 
for  shorter  pulses  and  it  induces  a  nonlinear  change  in  the  phase  slip,  which  is  known  to 
depend  on  pulse  energy  (pump  power). 

We  found  the  nonlinear  slip,  i.e.,  the  slip  induced  by  nonlinear  phase  and  nonlinear 
dispersion  effects,  to  be  well-approximated  by 

JNL  «  3 k"lc/T%s, 

where  k "  is  average-cavity  GVD.  tq  is  the  pulse  width,  lc  is  the  optical  cavity  length  and 
s  is  the  map  strength.  This  result  shows  that  the  phase  slip  that  is  induced  by  nonlinear- 
dispersion  vanishes  with  stronger  dispersion- management,  which  is  consistent  with  the  in¬ 
sensitivity  of  the  slip  to  pulse  energy  with  strong  dispersion-management .  An  additional 
effect  which  impacts  the  phase  slip  is  third-order  dispersion  (TOD),  which  can  be  mod¬ 
eled  bv  adding  [ikm / {^oP*T^)]uttt  (normalized)  to  the  right-hand  side  of  Eq.  (??  (with 
g  —  r  —  l  =  0),  where  km  is  the  average  TOD  coefficient,  i.e.,  the  net  TOD  per  round-trip 
of  the  cavity  and  P*  is  the  characteristic  power. 

Thus  the  laser  cavity  was  modeled  by  a  dispersion  and  nonlinearity  managed  nonlinear 
Schrodinger  equation  (perturbed  NLS),  that  takes  nonlinear  phase  (self-phase  modulation) 
and  small  nonlinear  dispersion  into  account.  In  our  work  wc  developed  a  detailed  asymptot  ic 
theory  which  gave  analytic  results  that  described  the  carrier-envelope  phase  shift  in  this 
dispersion-managed  NLS  equation.  Control  of  the  carrier  envelope  phase  shift  is  important 
in  applications.  It  is  a  key  aspect  to  obtain  highly  stable  optical  oscillators. 

Comparison  of  theory  with  experiments 

Motivated  by  the  question  of  whether  or  not  the  mode-locked  pulses  observed  in  the 
laboratory  were  actually  solitons,  a  number  of  experiments  were  performed  in  the  Dept,  of 
Physics  at  the  University  of  Colorado.  It  was  found  that  the  ‘pure”  dispersion-man  aged 
theory,  without  gain/loss  terms  (i.e.  in  this  case  g  =  r  =  I  —  0  in  Eq.  (??))  agrees  with 
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experiments  remarkably  well  Moreover,  in  recent  work  accounting  for  g,r,t  /  0  we  found 
that  when  the  system  mode  locks,  the  solitons  are  well  approximated  by  pure  dispersion- 
managed  solitons.  This  further  supports  the  comparisons  of  theory  and  experiment  and 
the  carrier  envelope  phase  slip  calculation  described  above.  In  Fig.  ??  the  full  width  half 
maximum  (FWHM)  is  plotted  vs.  the  average  pulse  energy  for  various  values  of  the  net  group 
delay  dispersion  (GDD).  The  dashed  lines  are  the  theoretical  values.  This  experimental 
research  also  shows  that  dispersion- management  soliton  concepts  are  broadly  applicable. 


Average  Pulse  Energy  ^(nJ) 


Figure  G:  Fundamental  pulse  parameters  in  a  mode-locked  Titanium: sapphire  laser.  The 
points  are  the  measured  temporal  FWHM  at  four  values  of  the  average  cavity  GDD.  The 
curves  are  the  solutions  of  the  Dispersion  Managed  NLS  equation.  The  legend  states  the 
GDD  values  in  fs2.  The  errors  for  the  GDD  are  approximately  1%  and  the  errors  for  the  r 
values  are  negligible  on  the  scale  shown  Inset  shows  the  breathing  dynamics  of  a  dispersion- 
managed  soliton  \u\2  as  it  propagates  along  z. 
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Noise  induced  linewidth  in  frequency  combs 

We  have  discussed  earlier  that  mode-locked  Ti:s  lasers  generat  e  trains  of  optical  pulses, 
whose  spectrum  consists  of  frequency  comb  lines.  These  combs  are  represented  by  evenly- 
spaced  frequencies  and  an  offset  frequency  that  is  proportional  to  the  carrier  envelope  phase 
(see  Fig.  ??). 

Random  physical  effects  can  induce  a  linewidth,  or  uncertainty,  in  a  comb  line  frequency. 
Some  random  effects  can  be  minimized  but,  as  far  back  as  1958,  Schawlow  and  Townes 
discovered  that  the  linewidth,  or  monochromaticity,  of  a  single-mode  continuous-wave  (cw) 
laser  is  fundamentally  limited  by  the  random  process  of  amplified  spontaneous  emission 
(ASE)  in  t  he  lasing  medium.  We  studied  the  limits  of  frequency  combs  associated  with  inode- 
locked  lasers  in  the  presence  of  ASE  noise  and  obtained  a  scaling  law  result.  In  particular,  we 
analyzed  the  frequency  combs  generated  by  trains  of  pulses  emitted  from  mode-locked  lasers 
when  the  center-time  and  phase  of  the  pulses  undergo  noise-induced  random  walk  which 
in  turn  broadens  the  comb  lines.  Detailed  asymptotic  analysis  and  computation  of  the 
ensemble-averaged  spectrum  revealed  a  time-frequency  duality.  In  particular,  t  he  increase 
of  the  standard  deviation  of  the  center-time  with  pulse  number,  and  the  increase  of  the 
linewidth  with  frequency,  are  closely  related.  More  precisely,  the  standard  deviation  of  the 
center- time  jitter  of  the  n’tli  pulse  is  found  to  scale  as  np^2 ,  where  p  is  a  jitter-exponent, 
while  the  linewidth  of  the  fc’th  comb  line  scales  as  k2^p.  The  linear-dispersionless  system 
(])  =  1)  and  pure  nonlinear  soliton  (p  =  3)  dynamics  in  lasers  are  found  as  special  cases  of 
this  time-frequency  duality  relation. 

Although  this  result  was  derived  for  mode-locked  lasers,  the  general  nature  of  the  result 
indicates  that  it  can  be  applied  to  the  stochastic  dynamics  associated  with  other  frequency 
combs. 
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Nonlinear  Optics  in  Waveguide  arrays  and  photonic  lattices 


Nonlinear  light  wave  propagation  in  photonic  lattices,  or  periodic  optical  waveguides, 
is  an  active  and  interesting  area  of  research.  This  is  due,  in  part,  to  the  realization  that 
photonic  lattiees  can  be  constructed  on  extremely  small  scales,  typically  a  few  microns  in 
size.  They  allow  the  possibility  of  manipulation  and  navigation  of  lightwaves  in  small  regions. 
Localized  nonlinear  optieal  pulses  which  oeeur  on  one  and  two  dimensional  backgrounds  have 
been  investigated.  These  backgrounds  ean  be  either  fabricated  mechanically  such  as  those 
comprised  of  AlGaAs  materials  or  all-optically  using  photo-refractive  materials  where  the 
photonic  structures  are  constructed  via  interference  of  two  or  more  plane  waves. 

In  two-dimensional  photonic  lattice  applications,  a  nonlinear  Schrodinger  equation  with 
an  external  potential,  derivable  from  Maxwell's  equations,  is  the  governing  equation.  The 
equation  in  normalized  form  is 

i*- — h  A u  —  V(x ,  y)u  +  \  u\22l  =  0  (2) 

(Jz 

Here  u(x,  y,  z)  corresponds  to  the  slowly- varying  complex  amplitude  of  the  electric  field  in 
the  plane  that  is  propagating  along  the  z  direction,  A u  =  uxx+uyy  corresponds  to  diffraction. 
I  r(x,  y)  is  an  external  optical  potential,  and  the  eubie  term  in  u  originates  with  the  nonlinear- 
Kerr  change  of  the  refractive  index  for  a  eubie  or  y®  nonlinear  index.  The  optical  potential 
corresponds  to  a  change  in  the  linear  refractive  index  of  the  medium,  which  can  be  achieved 
either  by  etching  or  with  an  all-optieal  induetanee  technique.  We  note  that  the  potential 
assumed  in  this  study  is  uniform  along  the  propagation  direction  2,  although  in  principle, 
non-uniform  potentials  in  2:  could  also  be  considered  in  a  manner  analogous  to  diffraction  or 
dispersion- management. 

Recently  we  have  begun  a  collaboration  with  faculty  in  the  Electrical  Engineering  De¬ 
partment  at  the  University  of  Colorado.  They  have  constructed  defeet,  dislocated  photonic 
lattice  systems  and  quasi-crystal  potentials  in  the  laboratory.  By  a  defect  we  mean  that,  only 
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one  or  two  lattice  sites  are  affected  by  eithei  removing  a  small  number  of  sites  or  enhancing 
the  amplitude  of  these  sites.  A  dislocated  photonic  lattice  system  has,  for  example,  a  line 
of  sites  merging  with  another.  A  quasi-crystal  structure  has  long  range  order  but  is  not 
periodic.  The  general  quasi-crystal  potential  we  study  is  given  by 

V(x,y)  =  ^ 

which  corresponds  to  the  diffraction  of  plane  waves  that  emanate  from  N  equally  spaced 
points  on  a  circle.  When  N  =  2,  3, 4, 6  regular  periodic  lattice  backgrounds  result.  The 
simplest  nontrivial  (i.e  nonperiodic)  case  is  N  =  5  which  is  often  referred  to  as  a  Penrose 
potential  this  is  named  after  R.  Penrose  who  discussed  such  functions  in  the  mathematical 
context  of  tiling  the  plane.  We  are  currently  studying  higher  order  and  more  general  cases. 
Defect,  dislocation  and  quasi-crystal  structures  are  often  seen  in  nature.  Typical  cast's  are 
given  in  Fig  ??  below. 


N- 1 

£• 


,2A:[cos(27rn/iV)a:+sin(27:n/iV)y 


vacancy  edge  dislocation 


y 


Penrose  (N=5) 


0*0  60^1 


X 


Figure  7:  vacancy  defect,  left,  dislocation,  middle,  quasi-crystal  order  N  =  5  Penrose 
structure-right 


In  our  research  we  constructed,  by  using  our  recently  developed  computational  methods, 
nonlinear  localized  pulses  which  are  associated  with  complex  defect  photonic  lattices  such 
as  the  ones  in  Fig  ??.  To  date  most  investigators  have  considered  regular  lattices.  Photonic 
lattice  systems  and  waveguide  arrays,  which  we  often  approximate  as  discrete  systems,  are 
natural  components  of  optical  design.  Since  they  can  be  fabricated  in  small  physical  envi- 
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ronments  they  allow  the  possibility  of  controlling  and  tuning  light  waves  in  confined  regions 
and  to  navigate  light  in  one  and  two  dimensional  networks. 

Dynamics  and  pulse  propagation  in  quadratic  nonlinear  optical  media 

In  many  applications  the  leading  nonlinear  polarization  effect  in  an  optical  material  is 
quadratic;  they  are  usually  referred  to  as  materials.  We  have  found  that  in  multidimen¬ 
sional  nonresonant  materials,  the  nonlinear  equation  governing  the  slowly  varying  enve¬ 
lope  of  quasi-monochromatic  wave  trains  is  a  coupled  nonlinear  system  involving  both  the 
optical  field  and  mean  terms.  They'  are  a  generalization  of  the  classical  nonlinear  Schrodinger 
equation.  We  call  these  equations  NLSM  systems  (M  stands  for  the  mean  contribution).  In 
water  waves  similar  scalar  systems  were  derived  in  1969  by  Benney  and  Roskes.  A  few  years 
later,  a  special  case  of  this  system  was  found  to  be  integrablc.  The  latter  system  is  frequently 
referred  to  as  the  Davey-Stewartson  (DS)  system. 

We  have  derived  both  scalar  and  a  vector  NLSM  systems  directly  from  Maxwell's  equa¬ 
tions.  The  vector  NLSM  systems  generalize  the  well  known  1  +  1  vector  NLS  equations  to 
multidimensions.  Such  vector  multidimensional  systems  are  new  in  mathematical  physics; 
they  have  not  yet  been  derived  in  other  physical  systems. 

In  our  research  we  found  that  the  y^2)  optical  system  can  exhibit  wave  collapse  corre¬ 
sponding  to  a  suitable  ranges  of  parameter  and  initial  data.  This  indicates  that  intense 
optical  pulses  can  occur  in  these  systems.  We  believe  that  experimentalists  will  be  able  to 
observe  this  phenomena  since  the  analagous  situation  was  recently  observed  in  cubic  non¬ 
linear  media.  These  theoretical  results  indicate  that  in  this  range  of  parameters  researchers 
must  be  careful  in  their  experiments  to  prevent  damaging  the  underlying  optical  crystal. 

The  NLSM  system  of  equations  possesses  nonlocal-nonlinear  coupling  between  a  field  that 
is  associated  with  the  first  harmonic,  with  a  “cascaded”  effect  from  the  second  harmonic, 
and  a  static  field  that  is  associated  with  the  mean  term  (i.e.,  the  zero  th  harmonic).  The 
system  that  we  analyzed  in  detail  can  be  written  in  the  following  non-dimensional  form. 
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(3h) 


(31.) 


where  u(i\y,z)  corresponds  to  the  field  associated  with  the  first-harmonic,  cor¬ 

responds  to  the  mean  lield,  and  v  and  p  are  real  constants  that  depend  on  the  physical 
paramet  ers. 

It  can  be  proven  that  the  above  system  can  admit  wave  collapse  of  localized  waves  when 
v  >  0.  We  remark  that  the  initial  conditions  are  u(x,y.Q)  =  uq(x<  y),  y,0)  =  do(x.y). 

such  that  equation  (??)  is  satisfied  at  z  =  0,  i.e.,  =  (|^o|2)j:- 

We  also  note  that  our  numerical  solutions  indicate  that  as  wave  collapse  occurs,  the 
solution  Lends  to  the  steady  state  mode  found  from  the  above  NLSM  system  (??).  The 
steady  solution  is  obtained  by  assuming  a  solution  of  the  form  u(x,  y,z)  —  F(x,  y)c,Xz  and 
<f>{x,yyz)  =  G(x.y ),  where  F  and  G  are  real  functions  and  A  is  a  positive  real  number. 
Substituting  this  ansatz  into  equations  (??)  gives 


A  F  +  ^A  F  +  F3  -  pFGx  =  0  . 


(4a) 


GXX  +  VGyy  ~  (F~)X  . 


(41>) 


We  plan  to  investigate  whether  localized  optical  pulses  can  be  obtained  when  we  have  an 
underlying  optical  lattice.  We  are  also  investigating  the  possibility  of  light  navigation  and 
optical  switching. 

Dispersive  Shock  Waves 

Shock  waves  in  compressible  fluids  is  a  classically  important  field  in  applied  mathematics 
and  physics,  whose  origins  date  back  to  the  work  of  Riemann.  Such  shock  waves,  which  wo 
refer  to  as  classical  or  viscous  shock  waves  (VSW’s),  are  characterized  b\  a  localized  steep 
gradient  in  fluid  properties  across  the  shock  front.  Without  viscosity  one  has  a  mathematical 
discontinuity;  when  viscosity  is  added  to  the  equations,  the  discontinuity  is  “regularized" 
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and  the  solution  is  smooth.  An  equation  that  models  classical  shock  wave  phenomena  is  the 
Burgers  equation 

ut  +  uux  =  vuxx  (5) 

If  v  =  0,  we  have  the  inviscid  Burgers  equation  which  admits  wave  breaking.  When  the 
underlying  characteristics  cross  a  discontinuous  solution  i.e.  a  shock  wave,  is  introduced 
which  satisfies  the  Rankine-Hugoniot  jump  conditions  which,  in  turn,  determines  the  shock 
speed.  Analysis  of  Burgers  equation  shows  that  there  is  a  smooth  solution  given  by 


1  1  ,  r  1  1 

“  =  2  "  2  tanl,(fe<;r  “  2!>) 

which  tends  to  the  shock  solution  as  v  — >  0.  Thus  the  mathematical  discontinuity  is  regu¬ 
larized  when  viscosity  v  is  introduced. 

Another  type  of  shock  wave  is  a  so-called  dispersive  shock  wave  (DSW).  Early  obser¬ 
vations  of  DSW’s  were  ion-acoustic  waves  in  plasma  physics.  Subsequently,  Gurevich  and 
Pitaevskii  studied  the  small  dispersion  limit  of  the  Korteweg-deVri.es  (KdV)  equation.  They 
obtained  an  analytical  representation  of  a  DSW.  As  opposed  to  a  localized  shock  as  in  the 
viscous  problem,  the  description  of  a  DSW  is  one  with  a  sharp  front  with  an  expanding, 
rapidly  oscillating  rear  tail.  The  Korteweg-deVries  (KdV)  equation  with  small  dispersion  is 
given  by 


lit  ~j-  UUX  —  t  Uxxr  (b) 

where  0  <  t  <<  1  regularizes  the  discontinuity  that  otherwise  would  be  present.  The 
mathematical  technique  used  to  analyze  DSW’s  relies  on  wave  averaging,  often  referred 
to  as  Whit  ham  theory.  Whitham  theory  is  used  to  construct  equations  for  the  parameters 
associated  with  slowly  varying  wavetrains;  it  provides  an  analytical  basis  for  DSW  dynamics. 
For  KdV  the  Whitham  equations  can  be  transformed  into  Riemann-invariant  form,  whic  h 
can  be  analyzed  in  detail.  One  finds  that  there  are  two  speeds  associated  with  a  DSW  one  is 
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the  speed  associated  with  the  frontal  wave  which  is  a  soliton  (located  at  j\s  in  the  figure),  hi  d 
the  other  speed  corresponds  to  the  group  velocity  of  near  linear  trailing  waves  on  the  rear  end 
(depicted  bv  X[  in  the  figure)  of  the  DSW.  The  picture  and  details  are  quite  different  from 
viscous  shock  waves  which,  for  example,  occurs  in  Burgers  equation;  of.  the  leftmost  1  ig.  ?? 
which  depicts  a  typical  DSW  associated  with  the  KdV  equation  and  a  classical  or  viscous 
shock  wave  (located  at  xc  in  the  figure)  associated  with  the  Burgers  equation.  Interestingly, 
the  structure  of  the  KdV  DSW  is  strikingly  similar  to  the  original  plasma  observations. 


Figure  8:  Left  figure:  typical  DSW  satisfying  the  KdV  eq.  (??)  and  a  classical  shock  wave 
satisfying  Burgers  eq.  (??);  middle  figure:  “in-trap”  blast  wave;  right  figure:  “out-of-trap” 
blast  wave.  Numerical  simulations  are  given  below  those  of  the  corresponding  experiment . 

Recent  experiments  in  Bose-Einstein  condensates  (BEC)  and  nonlinear  optics  have  en¬ 
hanced  interest  in  DSW  s.  The  BEC  experiments,  originally  performed  in  the  Physics  De¬ 
partment  at  the  University  of  Colorado,  motivated  our  studies.  Shock  waves  emanating 
from  a  blast  or  explosion  are  well  known  in  the  study  of  viscous  shock  waves.  Interestingly 
a  similar  situation  occurs  in  the  context  of  DSW’s.  The  dispersive  blast  wave  experiment 
and  computational  results  are  shown  Fig.  ??  middle  and  right  figures.  Recent  experiments 
in  nonlinear  optics  carried  out  in  the  laboratory  of  J.  Fleischer  at  Princeton  University  have 
also  observed  similar  blast-  waves  and  other  interesting  DSWr  phenomena. 


X 
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The  governing  equations  we  studied  are  a  defocusing  NLS  equation  with  an  additional 
external  potential.  In  BEC  this  equation  is  usually  called  the  Gross-Pitaevski  equation.  Im¬ 
portantly.  a  similar  equation  occurs  in  nonlinear  optics.  The  equation  is  given  in  normalized 
form  as 

ie^t  +  - V2T  -  V (r,  2,  t)V  -  |T|2T  -  0  (7) 

where  l/(r,  z )  contains  details  describing  the  cylindrically  symmetric  potential  trap  and  ad¬ 
ditional  laser  terms,  T  is  the  wavefunction  and  e2  is  a  small  parameter  related  to  t  he  conden¬ 
sate.  Two  different  configurations  were  studied:  the  so-called  ‘In-trap”  and  "out-ol-trap” 
cases.  The  two  cases  correspond  to  the  turning  on  of  a  tightly  focused  laser  beam  when  the 
BEC  was  in  either  a  trapped  or  an  expansion  (the  latter  is  called  nontrapped)  configuration. 
The  laser  beam  creates  a  dispersive  blast  wave  in  the  radial  direction. 

Numerical  simulations  were  found  to  agree  extremely  well  with  the  experiments.  In  Fig. 
??  (middle  and  right  figures),  experimental  and  numerical  results  are  shown  for  the  two  BEC 
configurations:  the  numerical  image  is  a  contour  plot  of  f  |T| 2dz  (darker  is  less  dense).  Our 
initial  analytical  studies  were  carried  out  on  the  1-D  semi-classical  NLS  equation  without 
the  external  potential  since  the  experiments  indicated  that  after  some  propagation  time  the 
DSW’s  were  approximately  1-D  in  character  and  weakly  influenced  by  the  potential. 

As  indicated  earlier,  an  analytical  approach  to  DSW’s  uses  Whitham  averaging  theory. 
Whit  ham  analysis  associated  with  NLS  equations  in  various  contexts  have  been  invest  igated 
by  a  number  of  authors  and  in  order  to  effectively  compare  the  analysis  with  BEC  experiment 
and  simulations,  by  ourselves.  Whitham  averaging  over  the  first  four  conservation  laws  of 
NLS.  using  a  1-phase  traveling  wave  solution,  leads  to  suitable  parameters  satisfying  a  system 
of  hyperbolic  PDE’s  which  can  be  written  in  Riemann-invariant  form.  Solving  the  Rieimum- 
invariant  system  corresponding  to  step  initial  data,  yields  a  description  of  an  NLS  DSW. 
The  NLS  DSW  is  a  slowly  modulated  wave  train  which  varies  from  a  large  trailing  wave 
which  is  well  approximated  by  a  moving  dark  (gray)  soli  ton,  to  a  nearly  linear  wavetrain  at 
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the  front  moving  with  its  group  velocity;  like  KdV  the  NLS  DSW  has  two  speeds.  The  1-D 
NLS  theory  was  applied  to  both  of  the  multi-dimensional  blast  wave  cases.  The  analytical 
results  were  very  good;  however  due  to  radial  and  potential  effects  there  is  a  difference  in 
phase  and  to  a  lessor  degree  in  amplitude. 

While  interactions  of  viscous  shock  waves  are  well  known,  the  situation  associated  wit  h 
DSW!s  is  still  at  an  early  stage.  We  have  made  some  progress  in  our  research,  but  more 
work  still  needs  to  be  done  in  order  to  develop  a  broad  and  detailed  understanding.  We  are 
investigating  DSW  interactions  in  physically  interesting  systems  by  employing  both  Whitham 
methods  and  asymptotic  analysis  applied  to  the  solution  obtained  by  the  inverse  scattering 
transform. 

We  are  studying  the  effect  of  higher  dimensions  on  the  interaction  of  DSW  s.  In  this 
regard,  we  mention  that  in  recent  nonlinear  optics  experiments  carried  out  in  .1.  Fleischer's 
laboratory,  interacting  DSWs  were  observed.  For  example  see  Fig.  ??  which  depict  various 
DSW  interactions. 


Figure  9:  Interactions  of  DSW’s;  top  (b):  nearly  one  dimensional;  middle  (d):  interacting 
cylindrical  waves  bottom  (f):  cylindrical  DSW  interacting  with  a  one  dimensional  DSW 
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We  also  analyzed  another  interesting  problem  in  the  theory  of  DS\V‘s:  the  piston  dis¬ 
persive  shock  problem.  Here  a  dispersive  shock  propagates  with  a  constant  speed  into  a 
dispersive  fluid.  The  one  dimensional  semi-classical  NLS  Eq.  with  imposed  moving  bound¬ 
ary  values  is  the  governing  equation.  For  small  piston  speed  the  result  is  a  DSW  separating 
the  fluid  at  rest  with  the  piston.  At  large  enough  piston  speed  there  is  a  bifurcation  of  shock 
behavior  and  a  locally  periodic  wave  train  is  generated  between  the  DSW'  and  the  piston. 

Dispersive  shock  waves  are  an  interesting  and  developing  area  of  research  which  we  believe 
will  play  an  increasingly  important  role  in  nonlinear  optics  applications  and  other  areas  of 
physics. 
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